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Abstract 

Scattering  of  a time-harmonic  electromagnetic  field  by  the  edge  of  a semi-infinite,  single- 
wall, zigzag  carbon  nanotube  (CN)  is  considered.  The  Wiener-Hopf  technique  is  applied 
to  determine  the  exact  solution  of  the  problem,  and  the  scattering  pattern  is  numerically 
calculated  in  the  vicinity  of  the  main  plasmon  resonance  frequency. 

1.  Introduction 

Since  the  discovery  by  Iijima  of  quasi-one-dimensional  cylindrical  crystalline  structures  of  car- 
bon atoms,  generally  referred  to  as  carbon  nanotubes  (CNs),  many  unique  physical  properties 
of  theirs  have  been  predicted  theoretically  and  detected  experimentally  [1].  In  particular,  with 
reference  to  their  optical  properties,  thin  films  comprising  aligned  CNs  have  been  described 
theoretically  as  composite  mediums  [3,  4]. 

A composite  medium  consists  of  a homogeneous  host  medium  with  periodically  or  randomly 
dispersed  inclusions.  The  inclusions  must  be  electrically  small  for  local  homogenization  to  be 
possible,  for  which  purpose  each  inclusion  is  represented  by  a polarizability  tensor  [5].  Fur- 
thermore, macroscopic  samples  of  a composite  medium  are  supposed  to  contain  a huge  number 
of  inclusions,  so  that  the  composite  medium  can  be  replaced  by  an  effectively  homogeneous 
medium. 

The  polarizability  tensor  of  a single  CN  in  isolation  has  been  treated  approximately  by  several 
researchers.  For  instance,  the  3-D  polarizability  tensor  of  a zigzag  CN  was  calculated  by  Ma  & 
Yang  [6]  when  its  length  L and  cross-sectional  radius  R are  small  compared  to  the  free-space 
wavelength  A = 2ir/k  (i.e.,  E « 1 and  kR  •C  1).  The  2-D  polarizability  tensor  (per  unit 
length)  of  infinitely  long  CNs  has  also  been  treated  [3,  4,  7].  However,  in  the  optical  frequency 
range,  the  typical  geometric  parameters  of  actual  CNs  satisfy  the  following  conditions: 

kR  <£.  1 , L R,  kL  r*j  l . (1) 

Such  conditions  are  characteristic  of  wire  antennas  at  microwave  frequencies  [8].  A wire 
antenna  cannot  be  characterized  by  a polarizability  tensor,  because  the  contribution  of  the 
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high-order  multipoles  to  the  scattered  field  is  strong  due  to  the  last  condition  in  Eqs.  (1). 
Scattering  by  a long  wire  is  much  too  complicated  to  be  expressed  via  a dipole,  and  arrays  of 
many  long  wires  can  not  therefore  be  homogenized  in  the  same  way  as  arrays  of  electrically 
small  inclusions  can  be  [5].  The  essential  quantity  required  is  the  scattering  matrix  (or  its 
equivalent)  of  a single  wire  [9].  From  this  quantity,  the  scattering  pattern  of  a wire  array  can 
be  calculated  [10].  Analogously,  the  key  problem  for  the  optical  response  of  CN  arrays  in  the 
optical  regime  defined  by  Eqs.  (1)  is  the  calculation  of  the  scattering  pattern  of  an  isolated  CN 
of  finite  length.  Of  course,  care  must  be  exercised  because  CNs  can  not  be  necessarily  assumed 
as  perfect  conductors  — unlike  wire  antennas. 

The  effective  boundary  conditions  for  a CN  are  non-trivial  [11].  Accordingly,  the  responses 
of  CNs  are  different  from  those  of  wire  antennas.  For  example,  strongly  attenuated  surface 
polaritons  [11]  and  plasmons  [2]  appear  in  CNs,  instead  of  the  weakly  attenuated  waves  of 
longitudinal  current  in  wire  antennas.  Yet  the  universality  of  macroscopic  electrodynamics 
means  that  certain  common  effects  are  possible.  In  particular,  we  expect  resonance  effects, 
which  can  arise  as  a result  of  the  interactions  between  the  edges  of  a CN. 

This  paper  addresses  the  electromagnetic  scattering  properties  of  CNs.  We  use  the  Wiener- 
Hopf  technique  [12]  and  the  effective  impedance  boundary  conditions  [11]  for  a semi-infinite, 
single-shell,  zigzag  CN.  The  scattering  amplitude  of  a finite-length  CN  can  be  expressed  in 
terms  of  the  scattering  response  of  a semi-infinite  CN,  using  the  edge-wave  method. 


2.  Theoretical  Framework 


Consider  a CN  aligned  parallel  to  the  z axis  of  a circular  cylindrical  coordinate  system  (r,  <p,  z) 
whose  origin  is  located  at  the  center  of  the  circular  cross-section  of  the  CN.  The  edge  of  the 
CN  can  be  either  closed  or  open.  If  closed,  the  edge  is  almost  hemispherical.  However,  the 
oxidation  of  CNs  makes  the  open-edge  configuration  more  probable  — which  is  fortuitous,  as 
that  configuration  is  the  more  suitable  of  the  two  for  theoretical  analysis.  But  the  scattered 
field  is  almost  independent  of  the  edge  configuration  if  the  first  of  Eqs.  (1)  holds  true,  in  direct 
analogy  with  hollow  and  dense  wire  antennas  [12].  Hence,  we  restrict  ourselves  to  the  open-edge 
configuration. 

Let  the  incident  field  be  E-polarized,  with  harmonic  time-dependence  of  e~iwt,  and  propa- 
gating at  an  angle  0q  with  respect  to  the  2 axis.  This  field  is  represented  by  the  Hertz  potential 
The  total  potential  V's  = where  corresponds  to  the  scattered  field  as  per 


„(s)  _ 1 d2tp^  1 d2ip^  1 ( d2ip^  ,2  .(sA 
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The  boundary  conditions  satisfied  are  as  follows  (in  Gaussian  units)  [11]: 
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(3) 


(4) 

Here,  $(<M)  = ^kazz^(R,4,  z)  sin2  0o  emerges  from  the  incident  field;  c is  the  speed  of 
light  in  vacuum;  the  parameter  f0  takes  spatial  dispersion  into  account  [11];  and  dzz  is  the  axial 
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conductivity  of  an  isolated  CN.  Both  dzz  and  Iq  have  been  calculated  via  quantum  transport 
theory  [11].  The  boundary  conditions  (4)  have  to  be  complemented  by  radiation  conditions  as 
well  as  by  edge  conditions  [12]. 

The  boundary  value  problem  may  be  solved  with  the  Wiener-Hopf  technique,  with  the  Jones 
method  employed  to  derive  Wiener-Hopf  functional  equation: 

J+(a)C(a)  = M«)  " 7^  • (5) 


Here,  $(<^>,  a)  is  the  spatial  1-D  Fourier  transform  of  <&((/>,  z)  with  a as  the  (complex-valued) 
spatial  frequency  corresponding  to  z;  Ki(-)  and  /;(•)  are  modified  Bessel  functions  of  order 
l > 0;  G(a)  = -r£-17~2,  £ = -4in azz/w,  T = 1 - l0a2k~ 2 and  7 = Va2  - fc2; 

while  J+(a)  and  'F-(a)  are  two  unknown  functions  to  be  determined  as  per  the  Wiener-Hopf 
technique.  In  order  to  solve  Eq.  (5)  analytically,  one  has  to  apply  the  usual  factorization  and 
decomposition  procedures,  and  the  exact  analytical  expression  for  the  scattered  Hertz  potential 
is  then  obtained  by  the  inverse  spatial  Fourier  transform. 

In  the  far  zone,  application  of  the  saddle  point  method  leads  to 


H 


(*) 


4S)  ~ Fi(e,0o) 


gifcVr^+z7 
fci/r2  4-  z'2  ’ 


(6) 


where  8 = 7r  — tan  l{r/z)  and 

■gf^fcgsin^)  cos  (0/2)  JijkR  sin  9)£^_ 

1 ’ 0 G+(fccos0o)(l  + cos0q)  X sin(0/2)  X (cos0  + cos0q)  G-(kcos0)  ’ 


(7) 


(•)  is  the  cylindrical  Hankel  function  and  ./;(•)  the  cylindrical  Bessel  function  of  order  /; 
while  G(a)  — G+(a)G-(a).  The  function  F}(0, 9q)  is  the  scattering  pattern  of  the  edge.  The 
full  scattering  pattern  of  a CN  includes  additional  components  accounting  for  surface  polaritons. 


3.  Numerical  Results  and  Discussion 

We  calculated  the  far-zone  scattered  power  density  P(9)  ~ \Fq(6,  0o)|2i  with  1 = 0 sufficing 
for  most  realistic  incident  fields.  Following  Ref.  [11],  we  set  the  inverse  relaxation  time  v = 
0.33  x 1012  s-1,  which  is  in  good  agreement  with  the  recent  measurements  of  dynamic  room- 
temperature  conductivity  of  mats  of  single-wall  CNs  [13]. 

Let  us  examine  P(6)  in  the  vicinity  of  /?  = 1,  where  (5  = tun/ 270  = 1,  h is  the  Planck 
constant,  and  70  is  the  so-called  overlap  integral  [11].  This  case  is  of  special  interest  since  it 
corresponds  to  the  main  plasmon  condition  [2]  for  all  types  of  CNs.  All  other  resonant  lines  are 
interpreted  as  its  satellites.  Sample  plots  of  P(8)  calculated  in  the  regime  0 < 9 < tt  — 0q  are 
presented  in  Figure  1.  Let  us  point  out  here  that  the  saddle  point  method  is  inapplicable  in  the 
vicinity  of  9 = 7r  — 6q. 

It  is  clear  from  Figure  1 that  a relatively  weak  deviation  from  the  exact  resonance  condition 
f3  = 1 leads  to  a significant  decrease  in  the  scattered  field  intensity;  and  we  deduce  thereby  that 
scattering  by  a semi-infinite  CN  is  essentially  due  to  plasmon  propagation.  We  also  infer  from 
Figure  1 that  forward  scattering  is  stronger  by  2 to  3 orders  in  magnitude  than  backscattering. 
This  effect  persists  for  all  types  of  CNs  in  a wide  frequency  range. 

To  conclude,  we  have  investigated  the  scattering  of  a time-harmonic  electromagnetic  field 
by  a semi-infinite,  single-shell,  zigzag  CN.  We  have  found  an  exact  analytical  solution  in  the 
framework  of  the  Wiener-Hopf  technique.  The  solution  found  will  serve  as  the  basis  of  the 
theory  of  light  scattering  by  single  CNs  and  CN-based  composites. 
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Figure  1:  Sample  plots  of  the  scattered  power  density  P(9 ) for  semi-infinite,  single-shell,  zizgag 

CNs  (m,  n = 0),  when  9q  = 7r/4.  For  comparison,  the  results  for  £ = oo  (i.e.,  perfect  conduction) 

are  also  shown. 
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